Equilibration of a one-dimensional quantum liquid 
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We review some of the recent results on equilibration of one-dimensional quantum liquids. The 
low-energy properties of these systems are described by the Luttinger liquid theory, in which the 
excitations are bosonic quasiparticles. At low temperatures the relaxation of the gas of excitations 
toward full equilibrium is exponentially slow. In electronic Luttinger liquids these relaxation pro- 
cesses involve backscattering of electrons and give rise to interesting corrections to the transport 
properties of one-dimensional conductors. We focus on the phenomenological theory of the equilibra- 
tion of a quantum liquid and obtain an expression for the relaxation rate in terms of the excitation 
spectrum. 

PACS numbers: 71.10.Pm 



I. INTRODUCTION 

The low energy properties of one-dimensional quantum 
systems are commonly described in the framework of the 
so-called Tomonaga-Luttinger liquid [ll-Q- This descrip- 
tion applies to systems of either bosons or fermions, but 
regardless of the statistics of the constituent particles, 
the excitations of the Luttinger liquid are bosons with 
linear spectrum. The Hamiltonian of a Luttinger liquid 
is given by 



nh 



ff - E «l9l^«^ + Jfi^NiN - Nof + vjJ^ 



(1) 



where hq is the annihilation operator of a bosonic exci- 
tation with momentum q propagating with velocity w, 
Q. The Hamiltonian ([T|) assumes that the system has 
a finite size L and periodic boundary conditions are im- 
posed. Apart from the occupation numbers of bosonic 
states, the energy of the system depends on two integer 
numbers, N and J . In the case of fermionic Luttinger 
liquids they can be interpreted in terms of the numbers 
of right- and left-moving particles as A^ = N^ + N-^ and 
J = N^' — N^ . Parameters ujv and vj have dimension 
of velocity and depend on the interactions between the 
particles; Ag is some reference number of particles in the 
system. The Luttinger liquid theory described by the 
Hamiltonian ([T|) has been successful in predicting a num- 
ber of interesting phenomena, such as the renormaliza- 
tion of impurity scattering in interacting one-dimensional 
electronsystems [13, [ill subsequently observed in exper- 
iments 12-15j. 

An interesting feature of the model ^ is the complete 
absence of coupling between the bosons. As a result, 
the lifetimes of bosonic excitations are infinite and the 
system does not relax toward thermal equilibrium. It is 
important to keep in mind, however, that Eq. ([1} is the 
exact Hamiltonian of the system only for the so-called 
Luttinger model 2] where the spectrum of the fermions 
consists of two linear branches, Cp — ±vpp. In a generic 
situation this is an approximation applicable only in the 
vicinity of the two Fermi points, and thus the Luttinger 



theory ([T]) applies only at low energies. In other words, 
Eq. Il]) represents a fixed point Hamiltonian in the renor- 
malization group sense, which should, in principle, be 
amended by additional contributions describing various 
irrelevant perturbations. The latter are the operators of 
third and higher powers in bq and 6j, which result in scat- 
tering of bosonic excitations. They adequately account 
for the curvature of the spectrum near the Fermi points, 
which gives rise to a multitude of interes ting phenomena 
studied in the last few years. (See Ref. [ig for a recent 
review.) 

Another important aspect of the Luttinger model is 
that the original fermions are classified as belonging to 
one of two species: the right- and left-moving particles. 
In realistic systems there is no fundamental difference 
between the particles moving in opposite directions, and 
a right-moving fermion may become a left-moving one 
upon scattering. These backscattering processes give rise 
to several interesting phenomena not captured by the 
Luttinger liquid theory. 

One example is the effect of backscattering on the 
transport properties of the system. Experimentally 
transport can be studied in quantum wire devices [17, 18], 
where a one-dimensional system is smoothly connected 
to two-dimensional leads. In the absence of interac- 
tions, the conductance of a quantum wire is quantized 
in units of e^//i, where e is the elementary charge and 
h is the Planck's constant. Interactions between elec- 
trons included into the Luttinger liquid theory do not 
affect conductance quantization (l9l - l21| . On the other 
hand, the backscattering processes excluded from model 
(dl) reduce the conductance [22]. More detailed theories 
of conductance of long uniform quantum wires relate the 
correction to conductance due to electron-electron inter- 
actions to the rate of equilibration of the electron liquid 

mill. 

The physics of equilibration of a liquid of one- 
dimensional fermions is the main subject of this paper. 
It is another example of a problem where backscattering 
processes are crucial. For particles with a realistic spec- 
trum, such as £p = p'^ /2m, the relaxation of the system 
to equilibrium involves backscattering processes chang- 



ing the numbers A^^ and N^ . On the other hand, in 
a Luttinger hquid the difference J = N^ — N^ is con- 
served, even if the irrelevant perturbations are taken into 
account. As a resuh, equihbrium states of the system 
described by the Hamiltonian ([ij are characterized by 
different chemical potentials of the two species of parti- 
cles, /i^ and /^^, and their relaxation to a single value /i 
is neglected. 

Below we discuss the mechanism of equilibration of 
one-dimensional quantum liquids beyond the Luttinger 
liquid approximation. An expression for the correspond- 
ing equilibration rate r~^ was obtained microscopically 
for the regimes of both weak [23] and strong [23| interac- 
tions. An alternative phenomenological approach [26ll27| 
based on the Luttinger liquid theory is applicable at any 
interaction strength and results in an expression for the 
equilibration rate t''^ in terms of the excitation spectrum 
of the system. The latter can be either measured exper- 
imentally or derived microscopically for specific models. 
In Sees. HIl - IVl we review the phenomenological approach 
|26l [27| and discuss the implications of the results for the 
equilibration rate to experiments with quantum wires. 



II. EQUILIBRIUM STATE OF A UNIFORM 
LUTTINGER LIQUID 

Let us first discuss the possible equilibrium states of a 
Luttinger liquid. In general, the equilibrium distribution 
is determined by the integrals of motion of the system 
|28l |. We will assume that the irrelevant perturbations 
resulting in weak scattering of bosons are added to the 
Hamiltonian ([l} . Then there are four integrals of motion: 
energy, momentum, and the numbers of right- and left- 
moving particles, iV^ = {N + J)/2 and N^ = {N-J)/2. 
The Gibbs probability of realization of a given many- 
particle state i is then given by 



^exp 



E, + uP, - ^i'^N^ - fi'-'N 



RatR 



(2) 



where Ei and Pi are the values of the momentum of the 
system in state i. To obtain the equilibrium distribution 
of Bose excitations one also needs the expression jSj for 
the momentum of the Luttinger liquid 



P^ppJ + J^qblbg, 



(3) 



where the Fermi momentum is defined via particle den- 
sity, pf = nhN/L. 

Using the expression 1^ one easily obtains the equilib- 
rium form of the occupation numbers of the boson states: 



iV„ = 



1 



.(v\q\-uq)/T _ l' 



(4) 



Note that as a result of momentum conservation the Bose 
distribution depends not only on temperature but also 



the parameter u, which can be thought of as the velocity 
of the gas of bosonic excitations. 

In addition to the bosonic occupation numbers, the 
state of the liquid depends on the zero modes N and J. 
According to Eq. ^ in thermal equilibrium the latter is 
peaked sharply near an average value J, which satisfies 



nh 



vjj 



UPF 



1 ^ 



(5) 



where A/i = /i — M • In a Luttinger liquid the ratio 
j = vjJ/L has the meaning of the particle current [8|. 
Expressing the latter in terms of the drift velocity Vd as 
j = {N/L)vd, we find 



Vd= u + 



A/i 
2^ 



(6) 



This expression shows that in an equilibrium state of the 
Luttinger liquid the gas of excitations moves at a velocity 
u different from the velocity Vd of the system as a whole. 
This decoupling is a result of conservation of J, which 
allows for the possibility of A/i ^0. In a realistic system 
the backscattering processes bring about relaxation of 
A/i to zero, and the velocities u and Vd equilibrate. 



III. EQUILIBRATION RATE 

In order to study the kinetics of equilibration of a 
Luttinger liquid one has to consider the corrections to 
the fixed point Hamiltonian ([T]). In the case of spinless 
Luttinger liquid the irrelevant perturbations are terms 
of third and higher orders in bosonic operators, such as 

^?i+g2^9i^92> b\^+q2-qJ>\iKi^'i2^ etc. Such perturbatious 
give rise to scattering of the bosonic excitations and to 
relaxation of their distribution function toward the equi- 
librium distribution (j4]). Since the scattering of bosons 
conserves their total momentum, the resulting distribu- 
tion is characterized by a velocity u, which can be easily 
obtained from the initial momentum of the whole gas 
of excitations. When the distribution approaches the 
equilibrium form ^ the typical scattering events involve 
bosons with energies of order temperature, and the scat- 
tering rate Tq scales as a power of T . For instance, in 
the case of a strongly interacting system the equilibration 
rate of the gas of excitations scales at Tq"^ ex T^ |29| . 

The backscattering processes required for the relax- 
ation of the velocity u toward Vd have been studied mi- 
croscopically in the regime of weak interaction in Ref. [22. 
Because of the constraints imposed by the conservation 
of momentum and energy the simplest non-trivial pro- 
cess involves three particles. Fig. [TJ The Fermi statis- 
tics requires that in the dominant backscattering process 
two particles are within the energy range of order tem- 
perature from the left and right Fermi points, whereas 
the third one is within T from the bottom of the band. 
As a result of such a scattering event the third particle 
backscatters, i.e., the numbers of right- and left-moving 



iiT- 




FIG. 1: In the model of weakly interacting fermions the domi- 
nant backscattering process involves three particles: one near 
the bottom of the band and the other two near the left and 
right Fermi points [23 ]. 



particles change by one. Since the backscattering parti- 
cle fills a hole deep below the Fermi level, the rate of such 
processes is exponentially small, r~^ ex e~-^^/^ |22l.[23l|. 
We will see below that the backscattering rate is exponen- 
tially suppressed at low temperatures for any interaction 
strength. 

The strong suppression of the backscattering rate 
means that at low temperatures the equilibration of the 
quantum liquid proceeds in two steps. First, the bosonic 
excitations come to thermal equilibrium with each other 
and their distribution function takes the form Q. This 
thcrmalization takes a relatively short time of order tq. 
Second, over much longer time r the backscattering pro- 
cesses equilibrate the zero mode J with the bosons. Dur- 
ing this time the velocity u of the gas of bosonic exci- 
tations approaches the velocity v^ of the liquid as the 
difference of the chemical potentials of the right and left 
movers A/j, relaxes to zero, see Eq. (JH). The time depen- 
dences of u and A/i should follow the usual relaxation 
law 



du u~ Vd 

'dt " ~' 



d ^ An 

dt T 



(7) 



Expression ([7]) gives the formal definition of the equili- 
bration time T. 

In order to study the relaxation rate r~^ at arbitrary 
interaction strength it is tempting to use the Luttinger 
liquid description of the system. However, this approach 
is incapable of describing the particles near the bottom of 
the band, Fig.[TJ which are crucial for the equilibration of 
the system. More precisely, the bosonic Hamiltonian ([1]) 
provides an adequate description of the excitation spec- 
trum of a quantum liquid only at low energies, namely 
|e| < D, where the bandwidth D ^ vpp- Indeed, for such 
excitations the spectrum can be linearized and consists of 
two independent branches, as required in the Luttinger 
model. On the other hand, any excitation with energy 
|£| ~ vpF is not accounted for by the Hamiltonian ([!]). 

This difficulty can be overcome as follows [26|. Be- 
cause the small probability of an empty state near the 
bottom of the band plays the crucial role in the physics 
of equilibration, let us first consider the spectrum of the 
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FIG. 2: Spectrum of a hole excitation in a quantum liquid. 
The states with energies below D are treated as excitations 
of the Luttinger liquid, whereas the higher-energy states are 
modeled as a mobile impurity. The hole can change its mo- 
mentum by 5Q = qi — §2 by absorbing a boson with momen- 
tum gi and emitting one with momentum q2. 



hole excitations. For non-interacting fermions a hole with 
momentum Q can be defined as an excitation of the sys- 
tem obtained by moving a fermion from state pp — Q to 
Pf- For a system with concave spectrum, such as the 
one in Fig. [U the hole represents the ground state of the 
system with the total momentum Q. We use this obser- 
vation to generalize the concept of a hole excitation to 
the case of arbitrary interaction strength, and define the 
hole as the ground state of the system with momentum 
Q. Because moving a fermion from one Fermi point to 
the other changes momentum by 2pF without changing 
the energy of the system, the energy Eq of the hole is a 
periodic function of momentum and vanishes at Q = 0, 
±2pF, ±4:PF, .... 

The holes with energies below the bandwidth D have 
a nearly linear spectrum. They are accounted for in the 
Hamiltonian ([TJ as superpositions of various bosonic ex- 
citations with the same momentum. The holes with ener- 
gies above D are not included in the Hamiltonian ([T]) and 
are treated as mobile impurities in the Luttinger liquid 
;30-38i. The exact value of the crossover energy scale 
D is not important as long as it is small compared to 
the maximum energy of the hole Sp^ ~ vpF and large 
compared to the temperature T. 

The mechanism of equilibration can be described as 
follows. For simplicity, we assume from now on that the 
liquid is at rest, Vd — 0. The gas of bosonic excitations 
equilibrates relatively quickly, and the occupation num- 
bers of bosonic states take the form ((4]), which applies in 
the region \q\ < D/v represented by two straight dashed 
lines in Fig. [21 In a generic case the total momentum of 
the excitations in the initial state was not zero, so the 
Bose distribution Q has a boost velocity u j^ 0. As a 
result of interactions between the bosons a small frac- 
tion of the particles are promoted above the energy D, 
where they are no longer described by the Hamiltonian 
(|T]). At arbitrary interaction strength the properties of 
these higher energy excitations are rather complicated, 
but the lowest energy excitation at a given momentum 
Q is the hole. Since eg ^ T, the occupation of the hole 



states is given by the Boltzmann factor 



where the probabihty current J has the form 



f{Q) 



- p-[^Q-nQ)/T 



(8) 



The presence of the correction ~uQ in the exponent is 
assured by the fact that the hole interacts and exchanges 
momentum with the thermahzed bosons. As a result of 
many such collisions the hole, with small probability, may 
increase its momentum Q above pp, after which it is more 
likely to fall toward Q — 2pF than return to the vicinity 
of Q — 0, see Fig. [5] As the hole approaches Q — 2pp, 
it enters the region of linear spectrum at eg < D, shown 
by dotted lines in Fig. [51 There it can again be viewed 
as a superposition of the bosonic excitations. 

As a result of this rare sequence of scattering events 
the bosons have transferred to the hole the momentum 
2pp. Due to the conservation of the total momentum 
([3]), this decrease of momentum of the gas of excitations 
means that the zero mode J — N^' — N'" has increased 
by 2, i.e., one fermion has been backscattered. Also, the 
decrease of the total momentum of the bosons means that 
the velocity u has also decreased, in accordance with the 
relaxation law ([7]). 

The equilibration proceeds very slowly because the 
hole must pass the point Q = pp in momentum space, 
where the occupation numbers are exponentially small. 
One therefore expects 



-i=Ce-^/^, A 



(9) 



To obtain the prefactor C, the kinetics of the scattering 
processes should be considered in more detail. 

We start by noting that the equilibration rate is con- 
trolled by a small region of momentum space near q — pp 
where the energy of the hole is close to the maximum, 
^ — ^Q ^T. The width of this region can be estimated 
as (m*T)-'^/^, where we introduced the effective mass of 
the hole as 



1 

m* 



(feq 



dQ^ 



(10) 



Q- 



Although the region is narrow compared to pp , it is wide 
compared to the typical change of the momentum of the 
hole in a single collision with bosonic excitations. In- 
deed, an elementary scattering event consists of the hole 
absorbing one boson and emitting another, see Fig. [5] 
Since the bosons are thermahzed, the typical change of 
Q is of order T/v, which is much smaller than {m*Ty^^ 
at T <C A. This estimate enables us to simplify the 
problem considerably. 

The motion of the hole in momentum space is random 
and occurs in steps that are small compared to the size 
of the critical region near the barrier. Such diffusion 
in momentum space is described by the Fokker-Planck 
equation [39j on the time-dependent distribution function 

fiQ,t), 



dtf 



dgJ, 



(11) 



J = - 



B{Q) 



f + 9q 



/. 



(12) 



Here prime denotes derivative with respect to Q, and 
B{Q) has the meaning of the diffusion constant in mo- 
mentum space. It is defined as 



B{Q) = Y}^Qf^Q,Q+SQ 

SQ 



(13) 



in terms of the rate Wq qi of hole scattering from state 
Q to Q'. 

In the steady state regime the Fokker-Planck equation 
(jlip is solved by demanding that the probability current 
J be independent of Q. To find the value of J one has to 
impose the boundary conditions on the occupation num- 
bers /(Q) on the two sides of the barrier. Considering 
that the size of the crossover region in momentum space 
{rn*TY''^ is small compared with pp, one can approxi- 
mate Eq. ([5]) as 



I{Q) 



-(eQ-upF)/T 



PF 



Q:^im*Ty/^. (14) 



This expression specifies the boundary condition upon 
f{Q) to the left of the barrier. To find one to the right of 
the barrier we notice that the hole states with momenta 
Q and Q + 2pp are identical, and the occupation of states 
with Q between pp and 2pp is given by Eq. ([5]) with 
Q^Q-2pp. This yields 

/(O) = e-(^«+"P-)/^, Q - PF » (m*r)i/2. (15) 



Solving the first-order differential equation (J12p with con- 
stant J one finds that the boundary conditions (|14p and 
(Uni) are satisfied for 



J = u 



PfB 



(27rm*T3)i/2 



-A/T 



(16) 



where we took the limit w — ?► and denoted B — B{2pp). 
A non- vanishing constant J means that {L/h)J holes 
are passing any given point in momentum space in unit 
time. Each hole moving from the vicinity of Q = 
to that oi Q — 2pp takes momentum 2pp out of the 
bosonic excitations. One therefore concludes that the 
total momentum of the bosons changes with time at 
the rate Pb — —2pp{L/h)J. Given that the momen- 
tum of the bosons distributed according to Eq. ^ is 
Pb — [ttLT^ /3hv^)u, we find ii — —u/t with 



1 



SB / vpp \ ' 



Pi 
2m*T 



1/2 



.-A/T 



(17) 



As expected, the equilibration rate has the exponential 
form ©. To fully evaluate the prefactor, however, one 
needs to study the hole scattering rate Wq^qi and obtain 
the diffusion constant (IT5)) . 



IV. HOLE SCATTERING RATE 

The scattering of a hole by bosonic excitations is a spe- 
cial case of the problem of dynamics of a mobile impurity 
in a Luttinger liquid [3l| . At low temperatures the lead- 
ing scattering process involves two bosons moving in the 
opposite directions. By absorbing one boson and emit- 
ting the other the impurity can scatter from state Q to a 
new state Q' without violating conservation of momen- 
tum and energy, Fig. [2] The authors of Ref. [3l| obtained 
the temperature dependence of the mobility of the impu- 
rity in a Luttinger liquid in this regime, /i oc T~*. Using 
the expression /i — T jB for the mobility (see [39], § 21) 
one concludes 



B = xr^ 



(18) 



The evaluation of the coefficient x presents an interest- 
ing problem. Microscopic calculations can be performed 
in the special cases of either weak or strong interactions 
|23l [25| . Interestingly, one can also obtain a phenomeno- 
logical expression for x iii terms of the spectrum of the 
mobile impurity (hole) in the Luttinger liquid [2 a . |27| . 
Here we review the latter approach. 

The diffusion constant B in the expression ([T7| for the 
equilibration rate should be evaluated at Q = pp. On the 
other hand, it is instructive to consider a more general 
problem and study the hole scattering rate Wq^q^^q in 
Eq. p^ for arbitrary Q. The latter can be found from 
the Fermi's Golden rule 

91,92 

y.8{eQ-£Qj^iQ^hv\qx\-hv\q2\), (19) 

where igi.gj is the matrix element of the process in which 
the hole absorbs the boson q\ and emits the boson (72, 
Fig. [2l Since the typical energies of the bosons are of 
order temperature, we will assume \8Q\ <^ pp- In this 
case one can easily obtain the momenta qi and 92 from 
the conservation laws, 



91 = ^^ 



gl^QI, .. 



1SQ + !^\SQ\. (20) 



Here Vq — Eq is the velocity of the hole with momentum 
Q. Using Eq. (|20l) one easily expresses the scattering rate 
as 



Wc 



,Q+SQ 



h?v 



^91(^92 + 1)1*91,92!' 



(21) 



To evaluate the matrix element ig^.^j we need to discuss 
the Hamiltonian of the Luttinger liquid in the presence 
of a mobile impurity. 

We start by writing the Hamiltonian ([Ij in an alterna- 
tive form M 



Ho = ^Jdx[KiV9f+K-\V<P)% 



(22) 



where the two bosonic fields 0(a;) and 6{x) satisfy the 
commutation relation 



[(t){x),V9{x')] =inS{x-x') 



(23) 



and the Luttinger liquid parameter K depends on the in- 
teractions between particles. The case of non-interacting 
fermions corresponds to K — 1. 

The Hamiltonian 1^^ can be brought to the form ([T]) 
with the help of the following expressions for the fields cj) 
and 9 in terms of the bosonic operators: 



V0(a:)= -^5]y^Msgn(g)(6, + 6t_^)e^^-/^ (24) 
9 

^^^'^^-'T.^l^i^-bUy'''^'- (25) 



The advantage of the form (1221) of the Hamiltonian is 
that the fields and 9 have clear meanings in terms of 
the observables characterizing the quantum liquid. For 
instance, the field 0(x) accounts for the fluctuations of 
the density of the liquid 



n{x) = hq -\ V4>{x) 



(26) 



where no — N/L is the average density [9] . Similarly, the 
field 9 is related to the momentum k of the liquid per 
particle 



k{x) = -nW9{x), 



(27) 



see, e.g., Ref. 1271. 

The coupling of the hole to bosonic excitations in the 
Luttinger liquid can now be obtained by considering the 
dependence £Q{n, n) of the energy of the hole on the den- 
sity and momentum of the liquid. Using Eqs. (|26p and 
p7| we expand £q(7t., k) in powers of the bosonic fields, 

SQ{n,K) = £Q{no,0) + -dneQ'V<j)-hdf,eQ'V9 

IT 

--dnd.eQV^V9 + ... (28) 

TT 

Here all derivatives of £Q{n, k) are taken at n = ng smd 
K = 0. Taking into account Eqs. (|24|) and (|25p one sees 
that the second-order terms in Eq. ([28l) contain contri- 
butions in which a boson qi is absorbed and boson 52 on 
the opposite branch is emitted. The respective matrix 
element has the form 



t. 



(a) 



\/kig2| 

2TThL 



9lr^Q' 



(29) 



where we assumed that the hole is at a; = and intro- 
duced the notation 



din 



Kdl 



(^ 
K 



dl 



(30) 



In addition to the terms coupling the hole to two 
bosons, Eq. (pS)) contains the contribution linear in 
bosonic operators: 



(31) 



where 



dL^VKd^-^d,, dn = VKdr, + ^d,. (32) 



The linear coupling terms ([31]) also contribute to the ma- 
trix element tq 
theory, 



but in the second order perturbation 



''91, «2 



/I 



91 '72 1 



2i:hL 



dLSq+qidREQ 



£q 



dREQ- 



^Q+qi 
q2dLSQ 



£Q~£Q-q2 -f|92| 



(33) 



where we assumed gi > and q2 < 0, i.e., positive SQ as 
in Fig. [21 It is important to account for the corrections 
to the momentum of the hole in the numerator which are 
due to the fact that the two perturbations of the form 
([3T|) act on the states of the hole with different values of 



The expression p3p can be simplified by taking advan- 
tage of the smallness of gi ~ (?2 ~ ^Q ^ Pf- The ex- 
pression in the brackets appears to scale as 1/SQ. This 
term is obtained by neglecting corrections to Q in the nu- 
merator and linearizing the denominators in gi and q2. 
However, for the specific values (l^n[) of the boson mo- 
menta the two contributions in the brackets cancel each 
other. Evaluating the next order terms in qi and q2 one 
finds 



/(b) - 



/I 



9192 1 



1 QlSq dR£Q 

2TrhL m*Q v + vq v — vq 



-dhVQ 



V~VQ 



Orvq 



dhSQ 



V + VQ 



(34) 



Here the momentum dependent effective mass of the hole 
is defined by l/mg = —Eq. The first term in Eq. ([34[) 
originates from the expansion of the denominators of 
Eq. (j33p to second order in qi and g2, whereas the re- 
maining two terms are obtained by accounting for linear 
corrections in the numerators. 

Finally, one more contribution to the scattering ma- 
trix element tq^^q^ is obtained when the hole couples to a 
single boson, Eq. pip, which in turn splits into two. The 
latter matrix element involves three bosons and should 
therefore originate from cubic in (p and 9 corrections in 
the Hamiltonian. For a fluid at rest the symmetry allows 
only for even powers of 9, so the correction must have 
the form 



Ha. = I dx[a0{y(j})iyef + a^{y(j)f 



(35) 
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FIG. 3: The three types of processes contributing to the scat- 
tering matrix element iqi.qa- (a) The first order process, in 
which the hole couples to two bosons, (b) The second or- 
der processes involving two perturbations, each coupling the 
hole to one of the two bosons, (c) The second order process 
where the hole couples to a virtual boson, which is separately 
coupled to bosons q\ and §2. 



The values of the coefficients ag and a^ can be related to 
the density dependences of the parameters v and K of the 
quadratic Hamiltonian (j22p by considering the correction 
to the total Hamiltonian Hq + Ha caused by a small 
change of particle density 5n. As a result one obtains 

M 



ae = —^dn{vK), 



a^ 



67r2^"(if)- 



(36) 



In order to find a contribution to tq-^^q^ we need the matrix 
element of Ha that absorbs a boson with momentum qi 
on one branch and emits a boson (72 on the other branch. 
Using Eq. ([55)) we obtain 



-l^=^r^^Vm<l2{qi - q2)\bl^bq^{bq,^q, +6;^_ ). 

(37) 
The second order calculation of the matrix element tg-^^q^ 
with perturbations ([37[) and pip yields 
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The three types of processes leading to hole scattering 
with absorption of boson qi and emission of boson 92 are 
illustrated in Fig. [3] Their total is given by 
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where 

Yq - dli^EQ 
-QrVq 



1 dL£Q dREQ 

m*Q V + VQ V - VQ 

OlEq vdnK ( dREQ 



dLVQ 



dREQ 



VQ 



VK \v- vq 



9lEq 

V + VQ 



VQ 

■ (40) 



An alternative way of evaluating the scattering matrix el- 
ement involves performing a unitary transformation that 
eliminates the linear coupling (l5lT) of the hole to the 
bosons [27| . Upon this transformation only the quadratic 



coupling remains, which is then evaluated in the first or- 
der, similarly to Eq. (f^ . The resulting expression given 
by Eqs. (49) and (50) of Ref.iSis equivalent to Eq. (p9l) . 
Using the expression (15^ for the scattering matrix el- 
ement in combination with Eqs. (|2T|) and ([T3)) one easily 
recovers the temperature dependence (jTS]) . The coeffi- 
cient X takes the form 



X 



4^ 



(41) 



Equations (fT7|) . (fT8|) . and (|4ip provide a complete ex- 
pression for the equilibration rate of a one-dimensional 
quantum liquid in terms of the spectrum of hole excita- 
tions and its dependences on the particle density n and 
the momentum per particle k. 



V. DISCUSSION 

In this paper we discussed the equilibration of a one- 
dimensional quantum liquid of interacting fermions. The 
conventional Luttinger liquid theory [8|, |9| of these sys- 
tems neglects the processes of backscattering. In many 
cases this is an excellent approximation since the corre- 
sponding scattering rates are exponentially small at low 
temperatures, Eq. (|17|) . However, the Luttinger liquid 
approximation does not enable one to treat a number of 
interesting phenomena in which the backscattering plays 
the crucial role. 

One example is the conductance of a long uniform 
quantum wire. The Luttinger liquid theory predicts per- 
fect conductance quantization in these devices, regardless 
of the interaction strength |19l - [2l| . On the other hand, 
it is easy to show that at weak electron-electron inter- 
actions a correction to conductance appears due to the 
backscattering processes [22, l23|. Interestingly, an ex- 
pression for the conductance of a quantum wire can be 
obtained for any interaction strength |24| . 



G 



h 



1 



rj.2 



3 v'^pi L + 2vT 



(42) 



The backscattering gives rise to a negative correction to 
the quantized conductance, which grows with tempera- 
ture and with the length of the wire L. In short wires 



the correction 5G 



is exponentially small, but it 



saturates at 5G ^ —{e'^/h){T/vpFY i^i long wires. 

Temperature dependent corrections to conductance of 
quantum wire devices have been observed in multiple ex- 
periments [4fl-[42| . The data shows excellent quantization 
of conductance at lowest temperatures and a negative 
correction developing as the temperature is raised. These 
observations are in qualitative agreement with Eq. (|42p. 
In comparing the data with theory it is important to 
keep in mind that our discussion so far has ignored spins, 
which appear to play an important role in experiments. 
The result ((42]) can be generalized to include spins |24| . 



but the evaluation of the equilibration rate of a system 
with spins is still an open problem. Another complica- 
tion is that most experiments study rather short wires 
which may not be treated as uniform. 

Our discussion of the equilibration rate did not as- 
sume Galilean invariance of the system. On the other 
hand, momentum conservation was assumed. Thus the 
results do not automatically apply to systems of interact- 
ing particles in periodic potentials, such as spin chains. 
In such systems umklapp scattering by the external po- 
tential may facilitate equilibration. On the other hand, 
electrons in GaAs quantum wires have an essentially 
quadratic spectrum Cp = p^ /2m, where m is the effective 
mass of electron in this material. Such electron system is 
Galilean invariant, which leads to a few simplifications. 
First, the Luttinger liquid parameter in this case is deter- 
mined by velocity of the bosons, K — Trhn/mv. Second, 
the dependence of the excitation energy on momentum 
K has the simple form [4j] 



eQ{n,K) =eQ{n) + Q—. 



(43) 



For momenta Q in the vicinity of pp one can expand 
£Q{n) = A(n) - (Q - ppf I2m* and find 



Y^,^K\ A" 



^A' 



(44) 



Upon substitution into Eq. (J4T|) one recovers the results 
of Ref. ^M for equilibration rate in Galilean invariant sys- 
tems obtained by a different technique. It is worth noting 
that in this case the equilibration rate is fully determined 
by the density dependences of the velocity of bosonic ex- 
citations V and the maximum energy of the hole A. 

Although our main focus was on interacting systems of 
fermions, the approach and the results should be equally 
applicable to systems of bosons. Similar techniques have 
been recently applied to dynamics of dark solitons and 
mobile impurities in bosonic fluids [37], |38| . Finally, it is 
worth mentioning that in integrable models apart from 
energy and momentum there are multiple additional con- 
served quantities, and one expects that no equilibration 
of the system should take place. In particular one should 
find r^^ — 0. This conjecture has been checke d pq for 
the Calogero-Sutherland [4J| and Lieb-Liniger [43 mod- 
els. More generally, one expects [23| that for integrable 
models the quantity Yq given by Eq. (|40)) should vanish 
for any Q. 
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